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EFFECT  OF  NONEQUILIBRIUM  PHONONS  ON 
SUPERCONDUCTING  STATES  WITH  TWO  COEXISTING  ENERGY  GAPS 

Andr/-Marie  Tremblay 

Laboratory  of  Atomic  and  Solid  State  Physics 
Cornell  University,  Ithaca,  NY  14853 

n  Gerd  Schon 

Institut  fur  Theorie  der  Kondensierten  Materie 
Universitat  Karlsruhe,  BRD 

ABSTRACT 

It  is  shown  that  under  the  influence  of  tunnel  currents,  super¬ 
conductors  may  exhibit  a  "first  order  transition11  to  a  state  with 
two  coexisting  energy  gaps.  An  energy  and  number  conserving  approxi¬ 
mation  for  the  collision  operator  is  used  to  explicitly  take  into 
account  in  the  theory  the  effect  of  nonequilibrium  phonons.  Quanti¬ 
tative  predictions  for  experiments  are  presented. 

1.  INTRODUCTION 


We  have  now  reached,  mainly  through  the  study  of  radiation 
stimulated  superconductivity  and  of  inhomogeneous  states,  the  stage 
where  many  nnonequilibrium  phases'*  (nonequilibrium  collective  ef¬ 
fects)  may  appear  in  superconductors,  immediately  bringing  to  the 
forefront  the  problems  of  stability  of  dissipative  states  which  are 
nowadays  under  intense  investigation!  in  the  field  of  nonequili¬ 
brium  statistical  mechanics  in  general.  Thanks  to  the  work  of 
Eckem,  Schmid,  Schmutz  and  Schon^  (ESSS)  we  have  now  a  large  num¬ 
ber  of  cases  where  the  stability  of  dissipative  states  has  been 
investigated,  including  some  examples  where  detailed  balance!  is 
violated. 

Although  the  work  of  ESSS  predicts  many  phenomena  which  can  be 
qualitatively  checked  with  experiment,  a  quantitative  comparison  is 
difficult,  mainly  because  phonons  were  assumed  to  remain  in  equi¬ 
librium,  a  condition  which  we  know  is  violated  in  most  experimental 
situations. ^ 


In  this  paper,  we  return  to  one  of  the  first  examples  (which 
is  now  part  of  the  more  general  picture  of  ESSS)  where  the  ideas 
of  Schmid^  on  the  stability  of  nonequilibrium  superconducting 
states  were  applied.  In  ref.  5  (hereafter  I)  it  was  shown  that  in 
superconducting  tunnel  junctions,  two  values  of  the  energy  gap 
could  simultaneously  be  a  solution  of  the  gap  equation.  The  re- 
suits  are  in  qualitative  accord  with  the  pioneering  experimental  j  f 

work  of  Dynes,  Narayanamurti  and  Garno6  and  of  Gray  and  Willemsen7. 

This  paper  reports  on  how  the  theory  of  I  is  modified  when  1 

nonequilibrium  phonons  are  included  in  the  simplest  possible  way.  /tkiUitj  °u 
This  allows  us  to  set  limits  on  the  domain  of  validity  of  the  theory  /' t 
of  I  and  to  make  quantitative  predictions  in  the  regime  where  we 
can  show  that  nonequilibrium  phonons  only  renormalize  the  tempera- 
ture  entering  the  Landau-Ginzburg  equation.  If  these  simple  theo -  /  ^ 
retical  predictions  can  be  checked  quantitatively  by  experiments,  / 
it  will  give  strong  support  to  the  theoretically  sound  global  1 
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picture  emerging  from  the  work  of  ESSS.  However,  if  a  quantitative 
comparison  with  experiment  is  not  successful,  the  theoretical  cal¬ 
culation  of  at  least  the  effect  of  nonequilibrium  phonons  will  need 
to  be  improved. 

2.  MICROSCOPIC  THEORY 

For  convenience,  we  first  summarize  the  results  already  obtained 
in  refs.  5  and  2.  We  consider  a  tunnel  junction  consisting  of  the 
superconductor  of  interest  (the  probe)  coupled  to  a  second  super¬ 
conductor  (the  injector).  We  momentarily  neglect  the  effect  of  non- 
equilibrium  phonons.  We  also  assume  that, 

a)  the  probe  is  thin  enough  that  we  may  assume  a  uniform  injection 
over  the  thickness  of  the  junction; 

b)  the  gap  of  the  injector  A*,  in  contrast  to  the  gap  of  the  probe  » 
is  not  appreciably  perturbed  by  the  tunneling  process.  This  can  be 
achieved  by  using  an  injector  which  is  much  thicker  than  the  probe. 
Qualitatively,  our  results  should  also  apply  to  a  more  symmetrical 
case?  but  since  it  is  hard  experimentally  to  make  the  injector  and 
the  probe  identical,  we  prefer  to  study  the  case  where  they  are 
markedly  different; 

c)  we  restrict  ourselves  to  temperatures  T  close  to  the  transition 
temperature  Tc  (Landau-Ginzburg  region),  (A,A^  <<  T) . 

The  equations  describing  such  a  system  have  been  derived  many 
times.  The  gap  equation  is, 

(a  -  BA2)A  «  -xA  -  S2V2A  (1) 

where  5 2  is  the  usual  Landau  Ginzburg  coherence  length,  a*(T  -T)/T  , 
(3)  /8tt2T2  and  the  "gap  control'*  (or  anomalous  term,  or  control 
function)  x  Is  defined  by, 

x  dE  |  ^(DtaCE)  (2) 

bOO 

where  6n(E)  is  the  angular  average  (isotropic  part)  of  the  energy 
dependent  deviation  of  the  distribution  function  from  its  local 
equilibrium  value.  Nj (E)  2  © ( | E | —A ) |e| /(E-A2)^  is  the  BCS  density 
of  states  and  £n(E)  obeys  a  Boltzmann  equation. ^  The  collision 
operator  Igp  can  be  written  down  in  the  relaxation  time  approxima¬ 
tion,  Iep[6n(E)]  -  Tg  6n(E)  where  can  be  taken  as  a  constant  tq 

equal  to  its  value  at  E  *  0  and  T  *  Tc  because  E  and  A  are  small 
with  respect  to  T.  To  order  (A/T)2  the  scattering-in  term  can  be 
neglected.  Finally,  the  term  representing  the  effect  of  tunnel  in¬ 
jection  can  be  derived  using,  for  example,  golden  rule  arguments. 
Since  only  the  part  of  6n(E)  which  is  odd  in  energy  contributes  to 
X  in  Eq.  (2),  we  can  write, 

(n£)D  -  2B  JV(E-eV)[nT(E-eV)  -  nT(E)  ]  +  (eV«-»-eV)j  (3) 

B  =  (8e2RftN(0))-1 

where  R  is  the  resistance  of  the  junction,  ft  the  (effective)  volume 
of  the  probe,  N(0)  the  normal  density  of  states  and  N*  the  injec¬ 
tor  BCS  density  of  states. 


It  was  shown  by  ESSS  that  the  most  stable  states  in  this  situ¬ 
ation  were  homogeneous  so  that  the  analysis  of  I  applies  under  more 
general  conditions.  The  global  stability  of  the  various  homogeneous 
stationary  states  can  be  determined  from  the  following  potential 
(generalized  free  energy)  which  can  be  obtained  from  the  stationary 
solutions  of  the  Fokker-Planck  eauation  obeved  bv  our  svstem. 


solutions  of  the  Fokker-Planck  equation  obeyed  b 
F  -  -2N(0) y'd3r j  dA'[a-PA'2  +  x (A * ) ]A * . 


obeyed  by  our  system, 
(A  * )  ]A f .  (4) 


Since  such  a  potential  function  (generalized  free  energy)  exists, 
our  system  obeys  the  principle  of  detailed  balance^-  and  its  be¬ 
havior  will  be  analogous  to  systems  in  thermodynamic  equilibrium; 
in  particular  it  will  exhibit  a  first  order  phase  transition  when 
two  minima  of  F  correspond  to  the  same  value  of  the  generalized 
free  energy  density®. 

To  proceed  any  further,  we  need  an  explicit  expression  for  x> 
the  gap  control.  The  Boltzmann  equation  for  the  distribution  func¬ 
tion  in  the  stationary  homogeneous  situation  reduces  to  6n(E)  * 

To^nE^D*  <luantity  X  is  comPuted  using  this  equation  and  Eqs.(2), 

eV 

(3).  The  first  order  expansion  in  ~zr  of  Eq.  (3)  is  plotted  in 

Fig.  1,  This  graph  of  fin^CE)  is  useful  to  understand  the  results 

given  in  analytical  form  in  I. 


Fig.  1:  The  first  order 
contribution  6n^1'(E)  to 
the  quasiparticle  dis¬ 
tribution  function  descri¬ 
bing  an  excess  number  of 
excitations  (particle¬ 
like  for  E  >  0  or  hole¬ 
like  for  E  <  0)  as  a 
function  of  energy. 


Consider  the  solid  line  in  Fig.  1  for  the  case  |eV|  ^  A^  +  A. 

We  can  qualitatively  understand  why  two  solutions  of  the  gap  equa¬ 
tion  may  simultaneously  be  possible,  as  follows.  Suppose  that 
A  *  A ,  is  a  self-consistent  solution  of  the  gap  Eq.  (1).  Then  the 
peak  in  6n(E)  at  I.IA^  will  be  convoluted  with  the  peak  in  N1(E)  at 
A  ■  A^^  and  will  give  a  large  negative  contribution  to  x  in  Eq.  (2). 
It  is1 thus  plausible  that  there  also  exists  a  solution  with  A>1.1A. 
since  in  such  a  case  the  peak  in  6n(E)  at  1.1A^  would  not  contribute 
to  x  a  more  positive  x  is  consistent  with  a  larger  gap  (see 

Eq.  (1)). 


The  effects  of  heating  on  6n(E)  will  in  most  instances  intro¬ 
duce  a  large  but  smooth  background.  If  spikes  similar  to  those  de¬ 
picted  in  Fig.  1  show  up  on  this  background  then  clearly  the  above 
qualitative  picture  of  why  there  may  exist  two  gap  solutions  will 
remain  valid.  On  the  other  hand  gap  enhancement  depends  very  much 
on  the  level  of  the  background. 

The  effect  of  nonequilibrium  phonons  becomes  important  when 
the  phenomenological  "escape"  time  Teg  becomes  much  longer  than  the 
phonon  scattering  time  xPh,  i.e.,  when  Teg  >>  tP^  where  xph  can  be 
expressed  as  a  function  of  tq  using  the  ratio  of  the  normal  elec¬ 
tron  and  longitudinal  phonon  specific  heats  C®1,  CP*1  respectively^, 


Tph/t0  =  SSeOKPb/CA/c*1), 


Using  the  energy  and  number  conserving  approximation  method  of 
Eckem  and  Schbny  we  find  that  the  effect  of  nonequilibrium  phonons 
can  be  represented  as  simple  heating  when 

1  «  T„/Tph  «  (T/a) (T/eV) , 

es  10 
In  such  a  case,  the  parameter  a  becomes  a1 

a'  *  a  +  6a,  =  a-2.59Bt  (eV/2T)2(T  /tph)  . 


(5) 


es 


(6) 


Note  that  6a,  does  not  have  any  sharp  step.  By  contrast,  x  computed 
with  has  a  step5.  It  is  basically  this  step  which  leads 

to  multiple  solutions  of  the  gap  equation. 

3.  SIMPLE  PREDICTIONS  FOR  EXPERIMENTS 
3.1  Qualitative  Aspects 

We  quote  a  few  qualitative  features  of  the  theory,  some  of 
which  have  already  been  experimentally  observed,  some  of  which  will 
hopefully  be  checked  in  future  experiments. 

a)  The  phase  transition  aspects  of  our  theory  have  already  been 
noted^»2  an(j  Seem  to  agree  with  experimental  findings^ *?:  at  a 
certain  voltage  (see  also  Sec.  3.2),  a  low  gap  region  (large  in¬ 
jection  current  density)  appears  in  the  probe  and  grows  relative  to 
the  larger  gap  region  (low  injection  current  density)  as  the  total 
injection  current  increases.  The  expected  hysteresis  has  also  been 
observed^* This  transition  in  the  junctions  is  analogous  to  a 
liquid  gas  transition  at  constant  pressure  where  the  relative  volumes 
are  controlled  by  the  total  volume.  The  total  injection  current 
which  controls  the  transition  has  no  analogue  in  any  of  the  other 
nonequilibrium  first  order  phase  transitions  studied  by  ESSS. 

b)  There  is  no  threshold  current,  conductance  or  quasiparticle  den¬ 
sity  in  our  theory.  The  two  coexisting  gaps  can  be  in  principle  ob¬ 
served  at  the  threshold  voltage  in  many  kinds  of  junctions,  even  at 
very  low  injection  current  density,  as  long  as  the  energy  differ¬ 
ence  between  the  two  coexisting  gaps  (see  following  section)  is 
large  enough  to  be  resolved  from  relaxation  time  smearing  effects 
and  other  similar  complications. 

c)  The  Chi  and  Clarke  gap  enhancement  (to  be  published)  and  the  co¬ 
existing  gaps  can  be  observed  in  principle  in  the  same  junction. 


d)  If  heating  effects  are  negligible ,  one  of  the  two  coexisting  gaps 
may  be  enhanced  with  respect  to  the  equilibrium  gap. 

e)  The  coexistence  of  a  superconducting  and  a  normal  phase  is  pos¬ 
sible,  as  can  be  seen  using  either  the  graphical  methods  of  I  and 
ESSS  or  an  analytical  method.  Since  this  possibility  is  probably 
harder  to  see  experimentally,  we  shall  not  give  any  more  quantita¬ 
tive  results  on  that  matter. 

Although  the  qualitative  features  described  in  paragraphs  a) 
and  b)  have  already  been  observed,  observation  of  the  phenomena 
described  in  paragraphs  c) ,  d) ,  e)  above  would  certainly  give  a 
stronger  experimental  support  to  our  theory. 


3.2  Quantitative  predictions 

Assume  that  Eq.  (5)  holds.  Then  Eq.  (6)  also  holds.  In  the 
homogeneous  case  and  when  |eV|  is  very  close  to  A  +  one  can  find 
a  simple  analytical  form  for  Eq.  (2).  The  transition  voltage  is 
determined  using  the  procedure  outlined  in  Eq.  (4). 

This  leads  us  to  our  first  quantitative  prediction. 


eV 


Ai  +  l/2(Ag<af)  +  A^(a')), 


(7) 


The  transition  voltage  eV  can  be  determined  experimentally.  It 
should  be  related  to  the  Injector  gap  and  to  the  experimentally  de¬ 
termined  values  of  the  two  coexisting  gaps  As(a*)  and  A.(af)  as  des¬ 
cribed  in  Eq.  (7).  Note  that  Eq.  (7)  and  Eq?  (8)  that  follows  are 
valid  only  when  the  condition  ^ 

!l o  /vrf)+M  h _ V  , 

a’  V  2T  /  Hi  («•)  /  1 

is  satisfied. 

The  temperature  and  conductance  dependence  of  the  two  coexis¬ 
ting  gaps  consitute  our  second  quantitative  prediction, 
roughly  constant  in  the  temperature  range  of  interest  and  Aj 
we  have 


As(a') 


If  A^  is 
'i  >> 


-1 


.-2 


Va’>  - 


A  (a ' )  'v* 


(Tc-T) 


3/4 


(Tc-T) 


7/4 


(8) 


-2. 


The  last  term  depends  on  the  square  of  the  conductance  (R  )  and  on 
T  and  by  assumption  is  smaller  than  the  term  proportional  to  the 
conductance. 


4.  CONCLUSION 

Five  qualitative  predictions  were  given  in  Sec.  3.1  and  two 
quantitative  predictions  in  Sec.  3.2.  Predictions  concerning  the 
absolute  values  of  the  observed  gaps  can  be  derived  but  they  are 
less  useful  because  they  contain  the  parameter  t  esA which  may 
be  harder  to  obtain  experimentally  Nevertheless,  it  is  useful  to 
observe  that  our  theory  is  not  inconsistent  with  the  fact  that  ex¬ 
perimentally,  the  two  coexisting  gaps  have  values  smaller  than  the 
equilibrium  gap.  If,  as  was  the  case  in  the  theory  presented  in  I, 
the  condition  Teg  /  Ag(a)+Ai’\  /  \>s 

~) 


is  realized,  then  the  largest  of  the  two  coexisting  gaps  has  a  value 
larger  than  the  equilibrium  value.  In  practice,  the  above  condition 
may  not  be  satisfied. 

It  would  be  interesting  to  extend  these  ideas  to  low  temperatures 
where  a  solution  of  the  deterministic  equation  also  leads  to  the 
possibility  of  a  phase  transition^  but  where  a  determination  of  the 
fluctuations,  which  we  must  know  to  determine  the  global  stability, 
is  more  complicated  than  close  to  Tc. 

If  the  predictions  presented  in  this  paper  are  verified,  we  are 
confident  that  our  theory  will  also  explain  new  experiments  which 
have  just  been  done^*. 

Finally,  guided  by  the  analogies  between  the  generalized  free 
energy  determining  the  stability  of  the  nonequilibrium  state  and  the 
usual  thermodynamic  free  energy,  one  can  speculate  that  future  ex¬ 
periments  may  exhibit  phenomena  analogous  to  spinodal  decomposition 
in  fluids^,  a  rapid  change  of  voltage  in  a  junction  is  analogous  to 
a  rapid  change  of  pressure  in  the  fluid.  The  current  and  junction 
conductance  are  analogous  respectively  to  the  volume  and  temperature 
of  the  fluid. 

It  is  a  pleasure  to  acknowledge  stimulating  discussions  with  V. 
Ambegaokar,  and  correspondence  with  K.  E.  Gray  and  B.  Huberman.  This 
work  was  supported  in  part  by  the  U.S.  Office  of  Naval  Research, 
Technical  Report  No.  3. 
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